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Abstract. We study the dynamics of a molecule's nuclear wave-function near an 
avoided crossing of two electronic energy levels, for one nuclear degree of freedom. 
We derive the general form of the Schrodinger equation in the n-th superadiabatic 
representation for all n € N, and give some partial results about the asymptotics 
for large n. Using these results, we obtain closed formulas for the time development 
of the component of the wave function in an initially unoccupied energy subspace, 
when a wave packet crosses the transition region. In the optimal superadiabatic 
representation, which we define, this component builds up monontonically. Finally, 
we give an explicit formula for the transition wave function away from the crossing, 
which is in excellent agreement with high precision numerical calculations. 
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1. Overview and main results 

We consider the dynamics of a nuclear wave packet as it travels through an avoided 
crossing of electronic energy levels. This problem has a long history in mathematics, physics 
and theoretical chemistry, and we feel that the reader should know both that history and 
some detail about the problem itself in order to appreciate our contribution. However, such 
a presentation would overload the introduction, and we shift it to a separate section. Here, 
we concentrate on giving a quick overview of our results, with little or no explanations. 
We will aim to give all necessary explanations later. 

We consider a time-dependent, two level Schrodinger equation of the form 

iedt7p{x,t) = Hi;{x,t), with H = -—d"^! + V{x). (1.1) 

Above, ip £ L-^(M, C^), / is the two-dimensional unit matrix, and the 2 x 2-matrix V{x) is 
the ('diabatic') electronic potential energy. We will usually write V in the form 

y(x) = p(x) f A. (1.2) 



(1.2) contains the implicit assumption that V is traceless. We assume that the system 
exhibits an avoided crossing of the adiabatic energy levels. Then a wave packet that is 
originally entirely in one adiabtic subspace undergoes non-adiabatic transitions when it 
travels through the avoided crossing region. In the adiabatic representation, these transi- 
tions are of order e globally, but they are exponentially small in e in the scattering region. 
It is known that there exist improved, superadiabatic, representations, such that the tran- 
sitions are exponentially small for all times; but no viable way of explicitly determining 



these was known. Our first main result is Theorem 4.4 , There we show how to compute, to 
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leading order in e, the n-th superadiabatic representation of the equation (1.1 ). The sym- 
bols of the off-diagonal coupling elements are obtained via a a recursive set of differential 
equations, given in Proposition |4.5[ 

The superadiabatic representations form an asymptotic series, i.e. they diverge for fixed 
e when n — > oo. Therefore it is natural to look for the optimal superadiabatic representa- 
tion, where the relevant quantities are minimal as functions of n for given e. To find it, and 
prove error bounds, one needs precise control on the large n asymptotics of superadiabatic 
coupling elements. We do not know how to obtain this control at present. Section |5] con- 
tains some partial results, including the asymptotics of superadiabatic coupling elements 



for large momenta, cf. Theorem 5.3 



In Section [6] we apply the previous results in order to compute non-adiabatic transitions. 
In contrast to the previous sections, the results of this final section are non-rigorous. The 
main result of Section [6] and in some respect of the whole paper, is an explicit formula for 
the exponentially small wave packet that makes the non-adiabatic transition in terms of 
data that are local in space and time. Let V'+,o (k) be the Fourier transform of the wave 
packet moving according to the Born-Oppenheimer approximation in the upper electronic 
surface at the time when its maximum reaches the transition point. If at that time one 
starts a wave packet with Fourier transform given by 

fJ{k) = sgn(fc)x{,2>45} sin (^) e-v I'^-WI (l + V^o'(^(A;)) (1.3) 

in the lower electronic surface and again evolves it according to the Born-Oppenheimer 
approximation, then a leading order approximation to the true time-evolution, i.e. to the 
solution of Schrodinger's equation (1.1), is achieved shortly after the wave packets leave 
the transition region. Here 

v{k) := sgn{k)\/ — 45 



is the momentum before the transition and 26 the energy gap. Formula (1.3) holds for the 
special case of constant eigenvalues ±5, i.e. p{x) = (5 in ( 1.2 ). It provides a practical way to 
correctly include non-adiabatic transitions at avoided crossings into the Born-Oppenheimer 
approximation. An analogous formula for Landau-Zener like avoided crossings will be the 
content of a future publication. 

While the derivation of (1.3) is not completely rigorous and many approximations enter, 
we show it to be in excellent agreement with high precision ab initio numerical calculations, 
for a wide range of parameters including moderately large e. Not only does it correctly 
describe the transition probability, but it yields the shape of the wave function after the 
transition to such a high accuracy th at, e.g. for the situation plotted in Figure |2l the 
relative error for ip- (k) from formula (1.3) and from a highly accurate numerical solution 

of the Schrodinger equation is of the order 10^^ uniformly in those k on which i/"- (k) 
is essentially supported. Put differently. Figure |2] shows both the true solution and the 
approximation by our formula (1.3). 

The method to derive (1.3) is to investigate the full time evolution of the transition 
wave packet in an optimal superadiabatic basis. In the time-adiabatic simplification of 
the problem, it is known [T, ^ that there exists an optimal superadiabatic representation 
in which transitions are (uniformly in time) exponentially small, and have the universal 
shape of an error function. For Born-Oppenheimer transitions, it is not immediately clear 
how to identify optimality of superadiabatic representations. In Definition [T] we propose 
a novel, natural criterion for optimality. We show that, if our criterion is met, transition 
wave packets in the Fourier representation have the universal error-function shape as 
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functions in time near their maximum. We also show that in a simple but important case 
our criterion can be met. 

The structure of our paper is as follows. In Section |2] we describe the basic ideas of 
the time-dependent Born-Oppenheimer approximation, and the main developments in the 
study of avoided crossings. In Section 3.1 we recall the construction of the adiabatic rep- 
resentation, which will help to understand the subsequent construction of superadiabatic 
representations in the remainder of Section [3] and Section |4j Section |5] is devoted to as- 
ymptotic results. Section [6] contains the applications to nonadiabatic transitions and has 
been discussed in detail in the previous paragraph. 



2. Introduction and history 

Time-dependent Born-Oppenheimer theory is the single most important tool for study- 
ing the quantum dynamics of molecules, just as the time-independent Born-Oppenheimer 
approximation |7] is for the theory of molecular bound states. The basic physical idea 
is that the electrons, being at least 2000 times lighter than the nuclei, move much more 
rapidly and thus quickly adjust their position with respect to the nuclei. In particular, if 
they start in the n-th bound state (for some fixed positions of nuclei), they should remain 
in the n-th bound state even though the nuclei are slowly moving; of course this will then 
be the n-th bound state with respect to the updated position of the nuclei. In turn, the 
nuclear quantum dynamical motion is determined by an effective potential given by the 
energy level of the n-th bound state of the electrons, as a function of nuclear position. All 
of this is expected to be true up to errors of order e, the ratio of masses between electrons 
and nuclei. 

Although named after Born and Oppenheimer, the time-dependent version of the theory 
was first proposed by Fritz London in ^23j . The mathematical investigation started with the 
work of Hagedorn [12], which made the ideas of London precise to leading order in e, under 
the assumption of smooth inter-particle potentials. Since then there has been considerable 
progress. In particular, the approximation has been pushed to arbitrary order in e, and 
extended to cover the case of an isolated subset of the electronic spectrum as opposed to 
just one single band j25l [28] [29] . 

Instead of trying to give a full account of all further contributions, we refer to the 
corresponding section of the excellent review article [18] by George Hagedorn and Alain 
Joye. 

The subject of this article is the time development of transitions between molecular 
energy levels, i.e. the dynamics of the part of the wave- function not obeying the Born- 
Oppenheimer approximation. These deviations from the adiabatic behaviour are of im- 
mense interest in quantum chemistry. In their most extreme form they occur when two 
energy levels of a molecule cross at a given configuration of nuclei, in which case Born- 
Oppenheimer theory breaks down completely and there are transitions of order one. This 
mechanism is now widely accepted as governing many important chemical reactions |30j . 
and is an active research topic both in the mathematical community |13| EI] and in theo- 
retical chemistry [U] fS] . Here, we will have nothing to say about it and instead concentrate 
on the related topic of avoided crossings of energy levels. The Born-Oppenheimer approx- 
imation then holds to leading order (and, as it will turn out, even beyond all orders in the 
scattering regime), but the remaining deviations are still of great interest. In fact, it is an 
avoided crossing situation that leads to the photo-dissociation of Nal, which is one of the 
paradigmatic chemical reactions in photochemistry [261 132] • 
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Let US consider a diatomic molecule, such as Nal. After discarding centre of mass move- 
ment and rotational degrees of freedom, of the nuclear degrees of freedom only the internu- 
clear distance x remains. For notational simplicity we ignore the spin degrees of freedom, 
and so the Hamiltonian is given by 

H = -^dl - ^A, + V^{x) + FenC^c, y) + V,{y): (2.1) 

where M is the reduced mass of the nuclei, ra^ is that of the electrons, y = (yi, . . . y„) are 
the positions of the electrons, and Ve and V^n are effective nuclear repulsion, electronic 
repulsion and nuclear-electronic attraction, respectively. We simplify (2.1 ) further by using 
atomic units (rrie = ^ = 1), putting M = e~^, and subsuming everything except the nuclear 
kinetic energy into the "electronic hamiltonian" Hc{x), acting as an operator in L-^(dy) 
for each x. As a result, ( |2.1[ ) reads 

H = -jdl + H,{x), (2.2) 

and on the time scale where a nontrivial nuclear motion can be observed the time- 
dependent Schrodinger equation is given by 

iedtip{x,y,t) = H'il;{x,y,t), ^l;{x,y,0) = i^o{x,y). (2.3) 

Let us now assume that the n-th eigenvalue of Hc{x), say En{x), is separated from the 
rest of the spectrum by a finite gap for all x, and non-degenerate. Let y i— > Xn{x,y) be 
the corresponding eigenvector in L^(dy) for each x. Then by Born-Oppenheimer theory we 
know that if we start the Schrodinger evolution with tpo{x, y) = (po{x)xn{x, y), the solution 
of ( |2.3| ) at time t will be given (up to errors of order e) by ip{x, y, t) = 0(x, t)xn{x, y), with 
(j){x, t) determined by 

iedt(t>{x,t) = -jdl^{x,t)+En{x)^{x,t), (t>{x,0) = (t>o{x)- (2.4) 

In this work we are not interested in the part that stays within the the energy band 
En, but rather want to study the orthogonal complement (in L^(dy)) of it, i.e. the part 
that makes the transition. To do so, we introduce a significant simplification of the model: 
we assume that the n-th electronic energy level interacts with one and only one other 
electronic energy level, say the (n — l)-st. The Hamiltonian for these energy levels alone 
is then given by 

H = -'jdl + V(.) with V(x)=(j.W fW^J, (2.5) 

where we assume that X(x) and Z{x) are analytic in a strip containing the real axis, and 
that p{x) = y^X'^{x) + Z'^{x) > c > 0. The latter corresponds to En{x) being isolated 
from En-i{x) in the full electronic Hamiltonian. H is now an operator in L^(dx, C^), and 
corresponds to the subspace of -L^(dy) spanned by the eigenvectors Xn and Xn+i- The 
approximation contained in (2.5 ) can be justified by arguing that if the remaining energy 
levels are separated from En and En+i by a gap that is uniformly larger than the minimal 
distance between En and En+i, then their interaction with the model system should be 
negligible. In general, we ought to be very careful with such claims, since the effects we are 
looking for are exponentially small. While it is conceivable that, by using space-adiabatic 
perturbation theory beyond all orders, one could rigorously justify that (2.1 ) gives the same 
transitions as (2.5) does in the relevant bands with exponential accuracy, there certainly 
exists no proof. Thus we think of (2.5 ) as an uncontrolled approximation and take it as the 
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starting point of our investigations. Note that even among the two-level systems, (2.5) is 
not the most general, since we forced the trace of V to vanish. While including a nonzero 
trace may lead to interesting effects, the case of vanishing trace is technically easier, and 
we defer the treatment of the general case to future work. 



It is instructive to simplify (2.5) even more by prescribing a ('classical') path x(t) for 
the nuclear motion of the two-level system, instead of considering its quantum evolution. 
The equation then reduces to 

iedt^it) = ( f ) m =: V{mt), (2.6) 

which is to be read as an equation for the occupation probabilities of the two relevant 
electronic energy bands given a certain motion of the nuclei. The limit e — > in equation 
( |2.6[ ) is known as the adiabatic limit. The adiabatic theorem states that any solution (j){t) 
starting at time to in an eigenstate ip{to) of the matrix V{to) remains an eigenstate of 
V{t) also at time t up to terms of order e as long as V{t) has two distinct eigenvalues. 
One says that transitions between the adiabatic subspaces are of order e in the presence 
of a spectral gap. However, in 1932, Zener |31j investigated an explicitly solvable instance 



of (2.6), namely X{t) = t/2 and Z{t) = 6/2. He observed that by solving (2.6) with 
an initial condition parallel to the eigenvector corresponding to +p{t) at t = — oo, the 
solution at t = -|-oo would have a component of magnitude e~'^^ Z^^^-* in the eigenspace 
corresponding to —p{t). In other words, in the scattering regime, the transition amplitude 
is exponentially small, much smaller than the bound obtained in the standard adiabatic 
theorem. Shortly after. Landau [20 argued that the same exponentially small expression 
should describe the scattering regime also for general analytic X{t) and Z{t) such that, 
at the minimum to of P^(0 = ^'^(0 + ■^^(^)) ^ is to first order approximated by the 
one considered by Zener. This gave rise to the famous Landau-Zener formula, which in 
itself has attracted much research, including [10] where more general situations leading 
to different prefactors where considered in generality for the first time, and [19] where a 
rigorous proof of the generalized Landau-Zener transitions was given for the first time. We 
refer to |17] for further information on the subject. 

Quantitatively the exponentially small scattering amplitude cannot be explained easily 
by any method involving just the adiabatic subspaces, i.e. the instantaneous eigenspaces 
of V{t). Indeed, when starting the evolution in one of the eigenspaces and monitoring the 
component 02 (t) parallel to the other one, one will typically observe a build-up of \4>2{t) \ up 
to order e, and later an eventual decay to the exponentially small final value. See |22| [6] 
for a numerical illustration of the phenomenon. One way to understand better what is 



actually going on is the complex WKB method: one solves (2.6) not on the real line, but 
on a curve in the complex plane that approaches the real line at ±00 and passes through 
the zeroes of the complex continuation of p{t). Then the solution is exponentially small all 
the way, and by this method Joye et. al. [19j prove the Landau-Zener formula. However, 
the method does not give any insight on what happens on the real line at finite times. 
This question was first investigated by M. Berry. In an influential paper [T] he expands the 



solution of (2.6) into a formal power series, and by truncating the resulting asymptotic 
series after n terms, he obtains a time-dependent basis of C^, called n-th superadiabatic 
basis. For all n, these bases agree with the adiabatic basis in the scattering regime, i.e. 
when the eigenspaces are approximately constant. Choosing n so that the remainder term 
in the asymptotic expansion is minimal. Berry shows that not only are transitions between 
the corresponding subspaces exponentially small, but they are also universal in the sense 
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that to leading order they are described by an error function for a wide class of matrices 
V{t). Later, Berry and Lim [2] refined these results by showing that some non-generic 
versions of V{t) would lead to different prefactors for the scattering amplitude. Both of 
these works were non-rigorous, and it was not until a decade later that Hagedorn and Joye 
jl5j succeeded to prove them in a special, non- generic case. Results in the same special 
case were obtained independently in [4 , but with a different method: Instead of expanding 
the solution of (2.6), the equation itself is transformed using adiabatic perturbation theory 



|29] . leading to superadiabatic representations, in analogy to superadiabatic bases. Each 
superadiabatic representation is given by a a unitary matrix Un{t) £ such that the off- 
diagonal elements of the matrix valued operator Un{t){iedt—V{t))U*{t) are of order e", but 
with a prefactor growing like a factorial. The optimal superadiabatic basis is the one where 
the off-diagonal elements are minimal, and there they are shown to be exponentially small 
Gaussian functions to leading order. Integrating the resulting effective equations using 
first order perturbation theory leads to Berry's results. It turned out that this method 
was sufficiently flexible to allow for a generalisation [5] covering all cases of interest, in 
particular the generic one and most of those appearing in [2 . So now (2.6) can be viewed 
as pretty much well understood. 

The situation is quite different for the true Born-Oppenheimer evolution (2.5). While 
it is known that it is possible (even for more general models) to construct exponentially 
accurate solutions |14| [25] [27] , all of these works merely give upper bounds on the error 
terms instead of computing them as functions of time to leading order, as was achieved in 
m for (2.6). Thus they say nothing about universality of transition histories, quantitative 
values of Landau- Zener scattering amplitudes, or scattering wave functions. Concerning 
the last point, there has been recent progress by Hagedorn and Joye [16]. They assume 
that the time evolution is started with a coherent semiclassical wavepacket, of sufficiently 
high momentum, located near x = — oo, and contained in one energy band. Under various 
assumptions, they prove that the portion of the wave-function making the transition to the 
other energy band is exponentially small and approximately Gaussian; they even give ex- 
plicit formulas for these wave packets in the scattering region, in terms of complex contour 
integrals. Of note, these formulas show that the exponential rate is larger than predicted 
by the corresponding Landau-Zener formula derived from (2.6), while the momentum of 
the transmitted wave-packet is larger than predicted by energy conservation. Intuitively 
this is due to the fact that fast portions of the wave packet are more likely to make the 
transition than slow ones, and the presence of those (as opposed to a perfectly sharp mo- 
mentum assumed in the approximation leading to (2.6)) increases both the exponential 
rate and the momentum of the transition wave. The method of their proof is again the 
complex WKB method. Since the results of [16] are of relevance to our asymptotic results. 



we will discuss the relation to our findings towards the end of Section 6.2 



3. Superadiabatic representations 



3.1. Representation as operators. Our starting points are (1.1 ) and (1.2 ). Switching to 
a superadiabatic representation means that we need to find a pseudo-differential operator 
Un in L^(M, C^) such that UnHU* is close to a diagonal (operator-valued) matrix in a 
suitable sense. As suggested by the notation, C/„ will be close to a unitary. The benefit of 
such a representation is that then, with = Un<p, the equation 



(3.1) 
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will decouple, up to error terms that we control, into two scalar equations. Moreover, 
an explicit control of the off-diagonal coupling terms in UnHU* will yield results on the 
superadiabatic transition histories simply through first order time-dependent perturbation 
theory. The most well-known instance of this procedure is the adiabatic transformation 
obtained by 

U _ ( cos(0(x)/2) sin(0(x)/2) \ 

^0^^^ ~ V sin(0(x)/2) -cos(0(x)/2) y • ^"^'^^ 
Since Uq diagonalizes V, we obtain 



e 



Since UqHUq acts on semiclassical wave functions that oscillate with frequency 1/e, the 
operator edx is actually of order one, and to leading order we get 



.2 



2^^"+ e'(x.) 



the standard expression for the adiabatic dynamics including the well-known derivative 
coupling. In particular, Hq gives the correct leading order dynamics inside the energy 
bands. However, if we are interested in exponentially small transitions, the adiabatic rep- 
resentation is not sufficient, and we need to find different representations (i.e. unitary 
transformations of L^(M, C^)) such that the powers of e in the off-diagonal increase. 

3.2. Symbolic representation. For this, we will be dealing with high order differential 
operators, and it is convenient to work in the symbolic representation. A full account of 
symbolic calculus can be found in |1H [M] . Here we only outline the main formulae that 
we are going to use, and do not touch the topic of symbol classes. Working in the symbolic 
representation means that we first replace x by g S M and ledx by an independent variable 
p E M in the definition ( 1.1 ) of i^. The factor e takes into account the semiclassical scaling. 
We then obtain 

2 

H{p,q) = \ + y{q). (3.4) 



where V is as in (1.2). The aim is now to find matrices Un{£,p,q) such that 

Hn{e,p,q) :=Un{e,p,q)#H{p,q)#Un{e,p,q) = diag + 0{e''+'). (3.5) 
Here # denotes the 

Moyal Product: For two symbols A{p,q) and B{p,q), their Moyal product is defined 
through A#B = with 

{A{p,q)#Bip,q))j = (2i)-^- J] ^ {d^d^A) [dp^s) . (3.6) 

The Moyal product is the natural product for semiclassical symbols and accounts for the 
non-commutativity of the operators they represent. It can be extended to symbols depend- 
ing on s, like Un, by representing C/„ as a formal power series Un{£,p, q) = "^j^ e^Un,k{p, q) 
and collecting powers of e, but we will only ever need the formula ( |3.6[ ) in explicit calcu- 
lations. 

Once Hn{£,p,q) is constructed, it can be recast as a pseudo-differential operator using 
the 



8 VOLKER BETZ, BENJAMIN D. GODDARD AND STEFAN TEUFEL 

Weyl- Quantisation: the operator W£{Hn) correspondmg to Hn{£,p,q) acts on a test 
function by 



1 

2TTe 



(3.7) 



With the connection between operators and symbols in place, we now need to actually 
construct the superadiabatic unitaries Un- For this we use the method of 
Superadiabatic Projections: that is, we seek symbols 



7r(")(e,p,, 



J2 ^^'^Ap^ i) (^ith TTj{p, q) e C 

j=0 



2x2^ 



such that 



F(p,g),vrW(p,g) 



# 



n+2\ 



(3.8) 

(3.9) 
(3.10) 



Above, [A,B]^ = A^B — B^A denotes the Moyal commutator. 

General theory ^9 guarantees the existence of vr^") for each n, and the same general 
theory states 



Lemma 3.1. There exists a semiclassical symbol Uni£,p,q) such that 

Uo{q)#Un{e,p,q) = Uo{q)Un{e,p,q) + 0{e) = l + 0(e), 

U:{e,p,q)#Un{e,p,q) = 1 + 0(e"+i) 

Un{e,p,q)#7r^''\e,p,q)#U:{e,p,q) = vr,, 
where tTj- is the projection onto the first component o/C^. 



(3.11) 
(3.12) 
(3.13) 



While the existence of the Un is known, they are usually tricky to calculate. Fortunately, 
there is no need to do so for getting the leading order of the superadiabatic Hamiltonian. 
This is the content of the following result. 



Proposition 3.2. Define Hn{£,p,q) as in (3.5), and assume the Un given in that formula 
fulfils Let Fn he given through Denote by 

ct{p,q) = {Uoiq)Fnip,q)Uo{q))i,2 

the upper off-diagonal element of UoFnUo, and let c~{p,q) be the lower off-diagonal ele- 
ment. Then 



Hn{e,P, q) 



pV2 + p{q) 



-ji+i „+ 



.n+l 



C„+l(P'9) 



+ 



0(e2) 0(e"+2) 
0(e"+2) 0(e2) 



Proof. The diagonal terms are immediate from first order adiabatic perturbation theory. 
For the off-diagonal terms, we multiply (3.13) with Un from the right and use (3.12) in 
order to find 

TTrUn = ?7„#7r(") #[7^ #[/„ = C/„#7r(") + e"+iC/„#7r(")#i? (3.14) 
for some symbol R, and similarly for UnT^v Iterating the above reasoning, we find 

C/n#vr(") = n,Un - e"+i?7„#7r(")#i? = n,Un#{l - e^+^R) + 0(e2"+2). (3.15) 
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We now use (3.14) in order to get 

(1 - 7T,)#Un#H#U:#7r, 

= (1 - vr,)#C/„#/7#7rW#C/:#(l + 0(e"+i)) = 

= (1 - 7rr)#C/„# (7rW#i/ + e^+'Fn+i + 0(e"+2)^ #^#{1 + 0{e^+')) = 

= (1 - ^r)#(^r#C/n + e''+'Un#7r^'^^ #R)#H#K#{1 + + 
+(1 - 7r,)#Un# + ^[7*^(1 + o(e"+i)). 

The next to last line above is ©(e^""*"^) by (3.15 ) and the fact (1— 7ri.)7ri. = 0, and multiplying 
with VTr from the right we find 

(1 - TTr)#Un#H#U*#7r, = s'^+^l - 7r,)Uo{Fn+l + 0{e))UoTT,. 

This is the result for the upper right off-diagonal element. For the lower left one, we 
interchange the roles of vTr and 1 — vTr, and obtain an additional minus sign from (3.10 ). □ 



Proposition 3.2 means that we can focus our attention entirely on superadiabatic pro- 
jections, and there will be no need to calculate superadiabatic unitaries. 



4. Superadiabatic projections 

Here we present the recursive scheme for calculating superadiabatic projections, getting 
more and more explicit as the section progresses. 

4.1. Matrix recursion for superadiabatic projections. It is easy to check that for 

we have ttqV = pttq. Hence ttq is the adiabatic (zeroth superadiabatic) projection, i.e. the 
projection on the eigenspace of V corresponding to +p. Projecting on the upper adiabatic 
subspace is just a matter of choice here, the same construction works when starting with 
the projection onto the other subspace. 

In order to obtain the higher superadiabatic representations, we note that 



Fn 



+1 



#,"+1 



[P /2 



fc=l 



Here, [A,B\^j^ = {A#B)k — {B#A)k is the coeffcient of [A,B]^ corresponding to e^. 
Using (3.6), we find 



n+l 



k=l ^ ' 



-\)\d\V) (5^7r„+i„fc) - (9,^7r„+i„fc) {d\V)\ . (4.1) 



Similarly, (3.13) gives 



n-l 



Gn+l = E 7rfc7r„+i_fc + ^ (vrfc#7r„_fc)^ + ^ {T^k#T^n-k- 1)2 + 

k=l k=0 k=0 

Finally, we can calculate Hn+i through 

Vr„+1 = Gn+l - TToGn - G„7ro -I- — TTo] . 

2p 



(4.2) 



(4.3) 
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The proof of (4.3) is given in |4|, Proposition 1, for the special case p = 1/2 and -F„+i = 
— i7r^(g). Since that proof apphes word by word to the general case, we do not repeat it 
here. Similarly, the important relations 

TToFnTTo = (1 - 7ro)F„(l - VTo) = (4.4) 

and 

(1 - 7ro)G„7ro = ^oCnCl - ^o) = (4.5) 

follow as in [3]. 

4.2. Transformed Pauli matrices. We need a more explicit recm'sive scheme in order 















I) 


' = ( 



to feasibly calculate the quantities F„, G„ and 7r„ given in (3.8 ) - (3.13 ). We now introduce 
such a scheme, following [3], but will use a different notation than [4] for reasons to be 
discussed below. Let 

;„-')■ -(j-i) («) 

be the Pauli matrices. We will need their representations in the adiabatic basis, given by 

'^x(g) = Uo{q)axUoiq), CTy{q) = Uo{q)ayUo{q), CTz{q) = Uo{q)azUo{q). (4.7) 

Above, note that Uq = Uq. The usual algebraic relations of the un-transformed Pauli 
matrices give 

fx'T'y = — ''■y'T'x = iCz, ''"x^'z = " '''zO'x = " i^y ) 

2 2 2 (4-8) 

O'yCTz = — CTzCTy = icTx, a^ = = tT^ = 1, 

where 1 is the unit matrix. Moreover, the special relations 

. (4-9) 

[c^xjTTo] = -io"y, [CyjTTo] = icTx, [cTzjTTo] = 0. 

can be easily checked. Here and henceforth, primes denote derivatives with respect to q. 

Remark: We will use the Pauli matrices as a basis to represent the 7r„ in. The basic 
idea is the same as in [4] , but there the basis matrices X, y, Z were chosen in an ad-hoc 
manner. It happens that the matrices X,Y,Z from |4] are linked to <Tx,cry,<Tz here, but 
unfortunately not in the most convenient way: We have X = icTy, Y = — ctz and Z = —tr^. 
This will lead to a serious clash of notation between the present paper and |1] , but we feel 
that the formulation of the problem in the widely used Pauli matrices justifies this. 

4.3. Superadiabatic projections through Pauh matrices. Let us first note that 
^z[<l) = y^o)/ p{<l)- This implies that the adiabatic projection is given by 

^o(p,g) = ^o('z) = ^(l + o-.((z)). (4.10) 

Indeed, vr*^ = ttq = ttq due io = Uq = 1, and 

[H{p,q),Mq)]# = Yil,n,{q)\ =eU,7,,{q) = '-^6' {q)a^{q). (4.11) 

L J # 1 / 

The final equality above follows from the fact cr'^{q) = —0'{q)crx{q)- In particular, we 
conclude that 

Fi{p,q) = l,pO'{q)cr^{q), 
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and Proposition |3.2| then gives 



Ho{e,p, q) 



pV2 + p{q) eip9'{q)/2 
-e\pe'{q)/2 p'/2-p{q) 



which is just the adiabatic representation (3.3) in symbohc language. 
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We now define the coefficients x„ through Wn by 

T^n{p,q) = Xn{p,q)(T^{q) + iyn{p, q)<Ty{q) + Zn{p,q)<r^{q) +Wn{p,q)l, (4.12) 

and emphasize that these coefficients have swapped names with those given in [4,. The 
prefactor i in front of will make for some more elegant formulas later, and slightly reduce 
the clash of notation with jj] . For our derivation of the recursions for the coefficients Xyi 
to Wni we ffist need to treat the derivatives of V appearing in (4.1): 



Lemma 4.1. We have 



= a„(g)cr,(g) + 6„(g)cr,(^), 



where an{q) and bn{q) are given by the recursions 

Mq) = p{q)^ boil) 
an+i{q) = a'^iq) + 0'{q)bn{q), bn+i{q) 





b'n{q)-0'{q)an{q). 



Proof. Use the fact V = ptr^ together with (4.9). 
The first step towards the coefficient recursion is 



(4.13) 
□ 



Proposition 4.2. Define Xn,yn, Zn and Wn through (4-12). Then 



n+l 



(2i)l! 



-b^ W< 



92 



(2i)22! 



a2yn-l + 



+ 



P 



(2i)33! 



+ ( |(z; + 0'x„) + 2 



(2i)l! 



92 



(2i)22! 



b2yn- 



q3 



(2i)33! 



+ 



V I 



d 

^^Ty^(ai2;n + biXn) + (^^p^(°3^n-2 + b'iXn-2) + 



1. 



Proof. This is just a calculation. The former parts of each bracket stem from the first term 
of (4.1 ), which computes to 

dqTTn = {x'n " O' Zn)cT^ + y^CTy + (^ + 0'Xn)(T^ + 



The latter terms of each bracket come from the various terms of second part of (4.1 ); for 
even k those are given by 



Qk 

Up 



(21)^=^! 



dqV,-Kn+l-k 



[akO-z + bkCrx,Xn+l-k(^x + il/n+l-k'^y + Zn+l-kf^z + Wn+l-k^] 



{2i)''k\ 



2i (iafe2/n+l-fc<''x + {bkZn+l-k - akXn+l-k)<^y - '^bkVn+l-kf^z 
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while for odd A; they are 



dp 

= ~ |-2i)fcfe; ^ {{o-kZn+l-k + bkXn+l-k)'^ + ak^n+l-kCTz + bkWn+l-kCT^) ■ 

Cohecting coefficients gives the result. 



□ 



It is remarkable that from this knowledge of -F„+i alone, we can derive a set of recursive 
differential equations that, together with zero boundary conditions at infinity, determine 
the coefficients x„ to Wn- 



Proposition 4.3. The coefficients Xn to Wn defined in (4-12) are determined by the fol- 
lowing recursive algebraic- differential equations: we have 

no) 



xi= zi= Wi=Q, yi = -ip 

4p{q) 

Moreover, 

yn = when n is even, Xn = Zn = Wn = when n is odd. 
For n odd we have 



Xn+l = - 

while for n even we have 

Vn+l = 



2p 



di 



2p 



1 / ^ (^'^ 



■J) \ ^ 











P 



di 



j=l V ) J- 



P 



(4.14) 
(4.15) 

(4.16) 

(4.17) 
(4.18) 
(4.19) 



Proof. For (4.14), note that 

1 



rp 1 'P^' ^ 1 P^' 



We now use (4.3) and Proposition 4.2 in order to prove the recursive formulae: from (4.5 ), 
we can deduce that [G„+i,7r„] = 0, and thus (4.9) implies that Gn+i is proportional to 
<Tz and 1. Consequently, the parts of 7r„-|_i proportional to cr^ and CTy arise from the last 
term of (4.3) alone, and comparison with (4.9) and Proposition 4.2 shows 

Xn+l 



1 P , 



dl 



2p \ i 



(2i)22! 



(2i)44! 



(4.20) 
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and 

Vn+l 



2p I " ' \ (2i)l! 



^ b^ w. 



iWn 



Q2 Q3 



(2i)22! 



(2i)33! 



71-2 



(4.21) 

Now, from (4.4) we can deduce that F„ is proportional to and CTy only, and so Propo- 
sition 4.2 immediately gives 



P 



(z; + e'xn) - 2 



(2i)l! 



92 



aiWn + 



p 



d„ 



(2i)22! 
53 



(2i)33 



,(4.22) 
(4.23) 



Now (4.15) follows inductively, and after that (4.16) - (4.19) are immediate from (4.20) - 
(4.23), noting that all we did is to add some terms that are zero, in order to get a more 
closed expression. □ 



We are now in the position to cast Proposition 3.2 into a more specific form, yielding 
our first main result. 



Theorem 4.4. Define Hn{e,p, q) as in Proposition 3.2, Xn{p, q) and yn{p, q) as in Propo- 

(4.24) 



sition 4-3, and put 



K^+l = -2p{yn+l±Xn+l). 



Then 



Hn{e,p,q) = yl + 



P 



n+l 



n+l 



-P 



+ 



0(e2) 0(e"+2) 
0(e"+2) 0(e2) 



Proof. From Proposition 3.2, we see that we need to calculate UoFn+iUo. By definition, 
Uoa^Uo = (Tx, and similarly for <Ty. Since is proportional to <Tx and cry, by comparing 



□ 



In order to make use of Theorem 14.41 we need to control the coefficients and Z'^. 
Since they are poynomials (of order n) in p, it makes sense to consider coefficients. We put 



Propositions 4.2 and |4.3| we obtain 

Uo{q)Fn+i{p,q)Uo{q) = 2p{q)yn+i{p,q)crx + 2ip{q)xn+i{p,q)cry. 



Comparison with (4.6) and Proposition 3.2 gives the result. 



Xn{p,q) = Y.P'""'<(l)^ 



(4.25) 



m=0 



with similar expressions for the other coefficients. Here, the index m in x'^ is an upper 
index rather than a power, and the choice of using for the highest power will turn out 
to be the most convenient one below. Differentiating now gives 



pXn+l-j 



n+l-2j 

E 

m=0 



(n + l-i-m)! ^„^i, 
(n + 1 — 2j — m)l 



2j~m m 



n+l— J ' 
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and thus the latter terms of expressions hke (4.16) are of the form 

L(n+1)/2J „+l-2j 



di 



(n + 1 — j — k)\aj 



n+l-2j~k k 



n+l-j 



n+l lm/2\ / I 1 I -M 



m=0 



Proposition 4.5. T/ie coefficients to w"^ defined in (4-25) are determined by the 
following recursive algebraic- differential equations: we have 



X™ = = u-i" = 0, m = 0,l, 



(4.27) 



Moreover, 



X. 



n+l 



1 / 1 

'2p 



(4.28) 



for n odd, while for n even we have 



Vn+l 



1 / 1 



n+l-m+j\( m-2j , m-2j n | 



[m/2j 



(4.29) 
(4.30) 

(4.31) 



Proof This simply uses (4.26) in (4.16)-(4.19). 



□ 



Remark: From the above equations, it is obvious that xj^ = = = w'^ = for odd 



m. But more is true. By induction and using (4.15), we find that 



_ m _ m 



w'^ = {) if m / 4A; + 2, 



(4.32) 



for k G No. 



5. ASYMPTOTICS 

It appears to be a very hard problem to determine the asymptotic behaviour of the 
quantities etc. as n — > oo, let alone prove it. So our results on this subject are rather 
incomplete, except in the case m = 0, i.e. for the term of the highest order (in p); in that 
case the asymptotics of x^(gf), y°(g) and z^{q) is known rigorously. 
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5.1. Highest order in p. For m = 0, the sums on the right hand side of (4.28)-(4.31 1 
are empty, and we retain 

4+1 = ^(2/^)', y^i = ^((4)'-^'^^), o = {z'j + e'yl 



After changing to the natural scale 



T{q)=2 / p(r)dr, (5.1) 







these are (apart from a change of notation) just the recursions appearing in the time- 
adiabatic case, which have been solved in [H |5]. 



The relevant results in that case read as follows: we introduce the natural scale (5.1 1 
and define f{T{q)) = f{q) for a given function /. Furthermore, we assume that 

^~e{T) = ^. '^ + e[{T), (5.2) 

dr T — ITc T + ITc 

for some 7 G M and Tc > 0, where 9^{t) has no singularities in {z G C : |Im(2;)| ^ Tc}, 
and only singularities of order smaller than one at itirc. As has first been observed in [2], 



(5.2) is the form of 9 for a large class of models, including the generic ones. For further 



discussion we refer to jH . 



From (4.27) we now conclude 



.9\q) . 9'{riq)) i ( ^7 



Let us write for the coefficient of belonging to p". Then by (4.24), we have 



^o,±(^^^ ^ / =F2p(g)x°(g) if n is even. 



-'^p{Q)yn{Q) if n is odd. 
Then by the results from ^ [5] , there exists /5 > such that 
4'±((?) = -ain)piq)i^{±ir+'dy{r) + in-l)\0{n-'')) (5.3) 

= a(n)p(g)i"(±i)"(n - 1)! ( . - . ) + (n - l)lOin~^) 

yl[q) = (n-l)!0(n-^). 

Above a{n) = '^'"^^J^^^ (1 + 0{l/n)) is the universal prefactor for the time-adiabatic tran- 
sitions. 

5.2. Estimates for lower orders in p. Unfortunately, the exact knowledge of the above 
asymptotics does not help us directly. The reason is that the terms p"'x^{q) and p^Unil) 
appear not to constitute the leading order contribution to Xn{p-, q) and yn{p-, g) as n — > cx). 
We expect this to be true in general, but it can be verified in the Landau-Zener case [3], 
where we have supg \x'^{q) \ ~ d^T{n+m/4) /T{n) for finite m and large n. We do not know 
the behaviour of when m is of the order of n; however, since numerical calculations 
clearly show that the latter case is where x"^ is maximal for fixed n, this is the regime 
that needs to be understood in order to do exponential asymptotics. 

In this work, all we can do is to prove a rough a-priori bound on the coefficients. While it 
is presumably not sharp, we will see that it identifies regimes of momenta for the incoming 
wave function such that the highest order in p determines the behaviour of the transition. 
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We recall the norms introduced in [4J, which we will use here. For Tc > and / C 
define 

-a+k 

d'^fit) ^ OO 



sup sup 



r(Q + A;) 

for a function / S C°° on the real line. We also define 

Fa,rSl) = {f^C^{I)--\\f\kl,a,r.)<Oo}. 

When Tc and / are fixed, we will simply write 

r(a + k) 



i, we 



(5.4) 



sup 

96/ 



a+fc 



J ll(/,",rc) 

and Ffj. In [1] we prove 

yk > 0, 



If II < 



/(r) dr 



< max 



(/,Q-l,rc) 



(a - l)|t - s\ 



(a) 



ll/9ll(/,a+/3,re) ^ B{a,(3) ||/||(7,„,^^) Ml (I ,p,r,) 



(5.5) 
(5.6) 
(5.7) 
(5.8) 



where B{a, /?) = r(a)r(/?)/r(Q + /3) is the Beta function. 

For the following a-priori estimate on the coefficients, we only treat the special case 
of constant eigenvalues, and without further loss of generality take p = 1/2. This is the 
case for which we will also derive the explicit transitions and transition histories later; the 
important Landau-Zener case will be treated elsewhere [3]. By (5.2), we have 6' G -^(i)- 

Proposition 5.1. Assume p = 1/2, and 9' G Fi{I). 

a) X™ G Fn{I) for any interval I, and the same holds for the other coefficients. 

b) Assume further that \I\ ^ min{rc/ ||0'||(;^^ , Tc/ ||0'||^j^^}. Then for every a > 1/2 
there exists Ca > such that for all m, n we have 



ll(/,n,Tc) ' II \\{I,n,Tc) ^ '^ct 



r(n + am) 



and 



"n ll(/,n,rc) ' II "^n \\{I,n,T^) ^ 



r(n) ' 

2Cq r(n + arn) 



1(1) 



r n 



Proof. We will proceed inductively and use Proposition 4.5, To this end, let us first note 
that by Lemma 4.1 and (|5.6l), (|5.8l), an easy induction shows that a„ and 6„ are in Fn, 



with ||an||(n) 5 ll^n||(„) ^ coln(n + 1) for some cq > 0. Thus, by Proposition 4.5 and (5.6) 



(5.8), we have Xn G Fn, and the same for all other coefficients. It remains to give bounc 
on the actual size of the norms. For m = 0, the bounds claimed in b) were given in H], 
cf. Theorem 2 there. Now let us assume that the claim b) holds up to m — 1 and up to n. 
Let us write 



Lm/2J 



1 



t (2i)^- 



n + l-m + j\ rn~2j 



m—2j \ 
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Then 



Lm/2J 



I cm II ^ 
Pn ll(n+l) ^ 



y - 



1 r(n + 2-m+j) T{n + l-j)T{j) 



(n + 2-m)r(j + l) r(n + l) 



X (116 



Lm/2J 



m— 2j 
Vn+l-j 



(n+l-i) 



+ a 



m-2j 
n+l-j 



(n+l-j) 



^ CaCo 1 + 



1(1), 



E 



1 In j r(n + 2 - m + j)r(n + 1 - j)T{n + 1 + am - (1 + 2a)j) 



2^^ J 



r(n + 2 - m)r(n + l)r(n + 1 - j) 



Clearly, the fraction of Gamma functions with j = 1 is the largest of all, and thus the 
factor 1/2-' allows to estimate the sum through twice its first term, giving 

,(n + 2 — m)r(n + am — 2a) 



ll(n+l) 



1) ^ CaC- 



where C does not depend on n or m. The same holds for all the other sums appearing in 
Proposition |4.5[ Using the recurrence relation there, we have 



|yn+l||(„+i) ^ Ikn \\{n) 



+ 



n 



1(1) II mil , ^ in + 2-m)T{n + am-2a) ^ 

II (n) + '-'a'^ ~ ^ ^ 



r n+1 



^ ,'T(n + am) T(n + am) ~(n + 2 — m)r(n + am — 2a 
^ C„ ( . — ^ + 2^- -r^ + C- ^ ^ 



a 



r(n) 

r(n + 1 + am) 



r(n + i) 

n + 2 



+ C 



r(n + i) 

(n + 2 — m)r(n + am — 2a) 



r(n + 1) \n + 1 + am ' r(n + 1 + am) 

The last term in the bracket above is 0{{n + am)~^"), and as a > 1/2, it vanishes faster 
than 1/n. Thus the bracket becomes smaller than one for large enough n. By choosing Ca 
so large that the induction hypothesis holds up to this n, we have shown the induction 
step for y^. The argument for is similar and simpler. As for z™, we have 



1 

||(„+i) ^ ^ 



,m/ 1 



(1) iFn ll(n) 



(n + 2 — m)r(n + am — 2a) 



and using (5.7), we see that 

1/ 



1(1) ll-^n ll(n) 



+ C 



r(n + i) 

(n + 2 — m)r(n + am — 2a) 



r(n + am) / ~ (n + 2 — m)r(n + am — 2a) 

^ L'a ~/ — ^ I 1 + C 



"n ll(n) 



r(n) \ r(n + am) 

The last bracket will be bounded by 2 for large enough n, and so the same reasoning as 
above shows the induction step for zJJ^. The proof for is similar and simpler. □ 



By piecing together intervals as given in Proposition |5 . 1 1 b) , and using (5.5), we obtain 

Corollary 5.2. For any compact interval I and any a > 1/2, there exists a constant C 
such that 

,r(n + am) 



sup|C(9)KC- 



q&I 



and the same for 
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5.3. Coupling function for high momenta. The final result of this section is about the 



asymptotic shape of the coupling functions given in Theorem |4.4[ for high momenta. 



From (4.25) it is clear that choosing p large enough will suffice to counter the growth of 



the coefficients as given in Corollary 5.2 To obtain precise statements, let us first note 
that by choosing 

Tr 



n 



ep 



we obtain 



eV4(g) = 2iV2e7rrcSin(7r7/2)e e ^-p-c e-p'^ {1 + OHepY^'^^")). 



(5.9) 



(5.10) 



This follows directly from the results in |4j . Recall also that under the assumptions 6' G Fi 



r^n^O ^, (n-l)! 

p for which this is also the leading order behaviour. 



and p = 1/2, we have p^x^ ~ p" ^"'^n'' at its maximum. We will now specify the regime of 



Proposition 5.3. Assume p = 1/2 and 9 E Fi, and / C M compact. Assume further that 
p = £^l^ , with 1/3 < /? < 1. Then there exists 5 > such that for all n ^ we have 



n , 



m=l 



Proof. By Corollary 5.2 and Stirlings formula, we have 



p"'{n 



m=l 



m=l 



'^{n)p" 



m=l 



for any a > 1/2. Clearly, this is largest for the maximal value n 
this case. Inserting e = p^^/^ into n = ^ gives 

J3_ 



pe 



SO it suffices to treat 



711-/3. 



Thus 



(n + am)""^ 

jinpm 



(n + am)"+"'^ 



n 



(3 13 
exp ( (n + am) ln(n + am) — (n H -m) ln(n) + m 



1-/3 



1-/3 



In Tr 



Now when /3 > 1/3, we can pick a > 1/2 such that /3/(l — (3) > a, and the exponent 
becomes negative for large enough n, and all m < n. The factor e"""* in the sum above 
then guarantees summability up to m = n without losing more than a constant, and the 
proof is finished. □ 



Together with (|5.10|), the previous result immediately gives 

We make the same assumptions as in 



Corollary 5.4. Define k 
Proposition 



0,± 



5.3 



and put n 



as in Theorem 
. Then 



pe 



f4i^{p,q) = ^2\yj2eTTTcSui{Ti-f/2)Q =p e 2.prc e=p^ (1 + C'(e'')), 
for some 5 > 0. 



(5.11) 
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This result is not completely satisfactory for several reasons. Firstly, in is only valid 
in the special case of constant eigenvalues p. However, an extension to general forms of 
the potential energy is merely a matter of, albeit laborious, routine. For the Landau- 
Zener model, this is addressed in a work in progress |3j. Secondly, our result only holds 
for momenta that scale with e, while ideally we would like to have asymptotics for fixed 
momentum. To get these appears to be a formidable problem, as knowledge about the true 
asymptotic behaviour of the for higher m is needed. We have put many efforts into 
that problem, with little tangible results. Thirdly, and most importantly from the applied 
point of view, it is not the asymptotic shape of the coupling functions we are interested in, 
but rather the shape of the time-dependent transmitted wave function. To obtain these, 
we will either need to solve the Wigner equation, or translate back into the language of 
opera tors using the Weyl quantisation. In both cases, it is inconvenient that the expression 
(5.11) contains terms of the form e"'' /^"^^^ , since these are in none of the usual symbol 
classes [23]. There is no theory, and worse, no calculus for dealing with such symbols. 

The way out of this dilemma is to not apply optimal truncation until the very end: we 
do a Weyl qunatisation for finite n, solve the corresponding PDE, and decide in the end 
how large we want n to be. This is the content of the final section. 



6. Transitions 



We now discuss transitions in the superadiabatic representations given in Theorem 4.4 
Since we do not have asymptotic information about the coefficients x^{q) for large n and 
m, we cannot treat the most desirable case of asymptotics in e for fixed momentum p. 



Therefore, except in Section 6.1 we will work in the spirit of the previous subsection and 
treat only large p. The approximation then consists in retaining only the terms in 
the more explicit formulas we provide. It will turn out that for those values of e where 
numerics of the exact Schrodinger evolution are feasible, our momentum does not have 
to be particularly large. These are also the cases where the result is not too small to 
be physically relevant. Concretely, we are talking about e between 1/10 and 1/50, and p 
between 2 and 5. We refrain from proving error estimates in this section, as they would 
be difficult to obtain and then could only hold for large p; rigorous treatment of such a 
restricted case would add something to our mathematical understanding, but not enough 
for us to consider it worth the effort. The value to practical applications of the formal 
calculations below is, on the other hand, potentially large. 



6.1. General transition histories. Our starting point is the Schrodinger equation in 
the n-th superadiabatic representation. Assuming that the potential V{x) approaches a 
constant matrix at spatial infinity quickly enough guarantees that all of the superadiabatic 
subspaces approach the adiabatic subspaces for large x. More precisely we assume for the 
moment that the limits lim^^-i-oo p{(l) exist, that 

lim = 

\q\^oo 

and that the limits are approached sufficiently fast. We now study transitions between the 
superadiabatic subspaces for solutions that are asymptotically, for t — > — cx), in the upper 
adiabatic subspace. More precisely, consider the Schrodinger equation 



\edtil^{x,t) = We{Hn{e,q,p))'tp{x,t) , 



(6.1) 
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where We{Hn) is the Weyl quantization of the symbol 







+ 



'71+ 1 







Here we split Hn{q,p), and in the same way also We{Hn), into a diagonal part of order 
one and an off diagonal "coupling" part of order e^~^^. Writing 

we can now use first order time-dependent perturbation theory in order to determine 
ip-{t,x) up to errors of order e^~^^: The perturbative solution of (6.1) with initial datum 
^p-^-{0,x) = 'tp-^-fl{x) and boundary condition limt^_oo ||V'-(t)|| = is 

V^+(t,x) = (e-i^+*V+,o)(^) + 0(^) 

and 

^P.{x,t) = -ie- /* (e-^(*-^)^- {WeK-^^) e-?^^+ V+,o) {x) ds + 0{e"+') . (6.2) 

J —oo 

Here 

i^n+iiP, q) = -2p(g)(-a;n+i(p, q) + Vn+iip, q)) 
is as given in Theorem |4.4[ and 



H±il>(x) 



^32 



In view of formula ( |6.2[ ), we need the Weyl-quantisation of n^^i{p, q)- Given that n^^i 
is a polynomial in p, we will be interested in a general formula for the Weyl quantisation 
of a symbol of the form p'^g{q). Moreover, we will later need the Fourier representation 
of the operator yVe(K~_|.;^). 

To this end we define the scaled Fourier transform 



1 



le \e 



1 



/27re 



--kx 



e £ 



f{x)dx. 



(6.3) 



Lemma 6.1. Let K{p,q) =p"''g{q) he a semiclassical symbol for some m € N. Then 



1 



/2tt£ 



g%k - ri) 



7] + k 



(6.4) 



Proof. We use (3.7) with 



V'(y) 



ip'^{r]) drj 



in order to get 

{WeK)i^{x) 



dyd^dr^g (^(x + y)) C e'M-y)+^y) 
dy di dr? g{y)C ei^^^+^^J^'^^^''-^^) ^'{r]) = 
dedryr(2(^-^))re^(2«-'')"^^(??) = 



(27re)3/2 
2 

(27re)3/2 
2 

1 

2^e 
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In the second line we changed variables y = (x + y)/2, and in the fourth ^ = 2^. We 
now apply the scaled Fourier transform to both sides of the above equations and use the 
formula 

^ ^' ^EbQ{i-r]-k) (6.5) 



dx e^(«-''"'=)^ 



in order to obtain the result. 

Remark: In position space, the Weyl quantisation of k(p 



□ 



(WeK)V'(x) 



(-i.rE(7)2 



= p"^g{q) is given by 
^{d^g){x){d"'-^i^){x), (6.6) 



which can be seen directly from (3.7) using integration by parts with respect to y, and 

(6:51). 



We will now give the momentum space version of (6.2) in a fairly explicit form. We write 
Q-s-H £qj, unitary propagator in the Fourier picture. We will also write k^^^ in the 
form 

L(n+1)/4J 
m=0 

according to our results from Section [4j Combining (6.2 ) and Lemma 6.1 now immediately 
shows that to leading order in e, and for any n G N, the component V'-.n of the solution 
to (2.5) in the n-th supeardiabatic representation is given by 



/27re 



Here, the operator Jn+i is given by Jn+i = Yl 



e^i^- J„+ie-^^^+^^,o)(^)ds 



L(n+1)/4J j4m^ ^-^^ 



m=0 



7] + k 



fiv)- 



(6.7) 



(6.8) 



6.2. Transitions for constant energy levels. Although general, (6.7) is not very help- 
ful when trying to calculate actual superadiabatic transition wave functions. To make 
further progress, we will make two simplifying assumptions. Firstly, we will treat the high 
momentum regime, as discussed above. Secondly, we will assume constant energy levels. 
This means that in (1.2) we will take 

p{q) = 6, and e'{q) 



17 



17 



(6.9) 



q-iqc q + iqc 

with 6 J- having no singular points of order greater or equal to one in the strip {z £ C : 
|Im (z)| ^ qc}. Constant p has the effect of trivializing at the same time the propagator in 
Fourier space. 



e 2e 



i;%k), 



and the transformation (5.1) to the natural time scale, with r(g) = 26q. Then, (5.2) 
becomes 

17 17 



dr9{r{q)) 



26q — iTc 25q + ir^ 



+ 0,{26q), 



with Tc = 25qc. 
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Now (4.24) together with (5.3) yield 

4'- ?a2a(n)(5(-i)"(n-l)!- 



1 



17 



17 



■(2(5)" V(9-i9c)" (g + igc)"/' 
Using the residual theorem, the scaled Fourier transformation in q of is given by 



(2(5)" 



2tt fkV^ ifci 
e E 



By (6.8), we find 



17 



A+if{k) = ^a{n + l) 



2\ " 



4(5 



{k + ri) e-^l'=-''l/(??) dr/. 



We now insert this into (6.7). For this it is useful to pick the origin of the time axis such 
that the transition occurs at t = 0. The evolution for ?/;_|_ is, to leading order, just the free 
Born-Oppenheimer evolution in the upper energy band; thus indeed we do not need to 
solve the full, coupled system to obtain the upper wave function at t = from the one 
given at some negative t and vice versa, which would defeat the purpose of what we are 
doing. With the above conventions, and our calculations so far, (6.7) reads 

vr:/(fc,o= ^"(^ + ^^ -i^(^v^'^)x 



ds / dr]{r] + k) 



k'^ — rf"- 
4(5 



e e 



(6.10) 



For the moment, we are interested in il)-^n {k, t) for t ^ 0. The integral in (6.10) converges 
as t — > 00, and in the limit can be calculated explicitly. To do this, consider the general 
integral 

" d??ds/(r/)efe(^'-'''-^^) = 
ds 



d7?/(?7)e^('='-'''-^^) + 



d??/(77)efe('='-'''-^'^) 



Setting fj = {k'^ — rj'^ — 4(5)/(2e) gives r] = ±(/c^ — 45 — 2fje)^^'^ and dr] = ^edfj{k'^ — 46 



.2 _ ^2 
2f]e)~^^'^, and so 



I 



+ = e 



d^ 



f{ + ^/W^^46^^2ffe) 
y/k'^ - 4(5 - 2fje 



df] 



j{ + ^/k^-45-2fje) 



^Jk'^ - 45 - 2f]e 



xjiv), 



where XJ is the characteristic function on J = {—00, ^"2^]- Similarly, 

f( - Jk^ -4S-2 



y'k^ - 4(5 - 2r]e 



We now use that J^2 dsdfjg(fi) e^^^ 
Fourier transform of g, to give 



drjdsf{r]) e^ 



^271^ J^dsg{s) = 27r(jf(0), where g is the inverse 



■^^^ fi+Vf^HfiW ^^ iffc2_45>0 
else. 
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Applying these calculations to (6.10) reveals that, for large positive t, we have 
TTja{n + 1) 



1p-,n {k,t) 



2 

+ ( 



-it(fcV2-<5) 



-1 + 



)e-Tl'=-v^^l^/(yp3^) 



(6.11) 



X{fc2>45}- 



The two terms in the square bracket in (|6.11 ) are clearly connected to positive and negative 

' ' £ 

incoming momenta, respectively. The second line will be negligible if either A; < or ifip^^o 
is concentrated on the negative half axis, while the third line will be negligible if /c > or 
if V'+,o is concentrated on the positive half axis. This shows the intuitively obvious fact 
that the transmitted wave packet will travel in the same direction as the incoming one. It 
also shows that we can replace \k it \/k'^ — 46\ with — Vk'^ — 46\ without changing the 
leading order result. We further streamline ( |6.11 ) by replacing a(n + 1) by its asymptotic 

value ^^^^[^j^, and by introducing 



v{k) = v{k, 6) = sgn(k)\/A;2 _ 4j_ 



Finally, we note that asymptotically, all the superadiabatic subspaces agree, so that (6.11) 
actually gives the asymptotic adiabatic transition. We thus conclude: 

After the transition, the wave function in the initially unoccupied adiabatic subspace is 
given by 



sgn(k) sin (J^ 



i^- {k,t) 

A few comments are in order 



it{/cV2-<5) ^-^\k~v{k)\ 



1 + 



flyl ^+,0 {v{k)). 



(6.12) 



1) The occurrence of the indicator function X{k^>i5} can be interpreted in terms of energy 
conservation: Any part of the wave packet that makes the transition obtains, in addition 
to its kinetic energy, the potential energy difference between the electronic energy levels. 
Thus, there cannot be any kinetic energy smaller than 26. Recall also that these 
transitions are radiation-less: instead of being radiated away from the molecule in the 
form of a photon, the energy is transferred into kinetic energy of the nuclei. 



2) We can read off a momentum shift from (6.11 ). We assume that V'+,o is a semiclassical 
wave function, and write iV'+ol^)! ~ e~^(^y'=^ . Let us assume for convenience that the 



Q-M{k)/e 

absolute minimum k^, of M{k) is on the positive real line. We find 

1 

l9c 



\n\i^n{k)\ « --{q^{^/v^ + A5-v)+M{v)). 



(6.13) 



Purely by energy conservation, one would expect the transition wave packet to be maximal 
when V = k^. However, since v ^ — V Is decreasing, the minimum of ^/v"^ + 45 - 

V — M{v) is shifted to the right. One can quantify this shift when M is given explicitly. 



3) The last point also shows that the term X{fc2>4<5} in (6.11) is of little consequence in 
practice. Since A;^ > 4(5 is equivalent to w > 0, and since only a small region around its 
maximum matters for the semiclassical wave function ij)- „ , we can safely leave out the 



factor X{A;2>4<5} in (6.11) without changing the leading order result. 

4) ( |6.11 ) depends on n only through the convergent prefactor a{n). In particular, we do 
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not need to know the value of the optimal n in order to obtain the correct leading order 



transitions. But in order to justify (6.11) it is of course important to choose n such that 



the error terms are smaller than the leading term. Moerover, the n-independence of the 
leading order term is a special feature of the constant eigenvalues, and cannot be expected 
to persist in more general models. 



5) Let us compute the transition rate from (6.12), in the limit of large momentum and 
small momentum uncertainty. We choose 

C 



tp+,0 (k) = exp 



2e 



{k-PoY 



When C is large, the minimum in (6.13) is taken very close to v = po, implying that 
k* ~ \/pq + 45. The value of the exponent at the maximum is then given by 



■7' 



p^ + 4:6- Po) 



which is the transition probability for momentum pQ. For large po, we have \/pq + 4:6— po 
25 /pq, so that the transition probability in this regime is given by 



exp 



2qcd 
Po£ 



exp 



PoeJ 



The latter is precisely the Landau-Zener transition probability for the parameters chosen 
in ( |6.9[ ), cf. e.g. [6], where one has to replace e by epo throughout. 

6) Another formula for the asymptotic shape of a non-adiabatic scattering wave function 
was given (and proved) by Hagedorn and Joye in |T6]. We do not give the formula here, as 
we would have to introduce too much additional notation, and refer instead to Theorem 
5.1 of [16]. But we comment on the differences between their result and ours. The first 



difference is that while (6.12 ) is for constant p, the work of Hagedorn and Joye covers the 



Landau-Zener situation. So, a direct comparison is not possible at present, but a version 



of (6.12) for Landau-Zener transitions is work in progress. The advantage of the result 



by Hagedorn and Joye it is that is rigorous, and that it covers the asymptotic region 
of arbitrarily small e. In contrast, by the arguments of Section |5] (6.12) is probably not 
asymptotically exact, although it works excellently for all cases that we were able to test 



on a computer; see below. The great practical advantage of (6.12) over the results of [TB! 



is that it contains no complex contour integrals, and that it does not rely on a second 



order approximation to the incoming wave function. Indeed, we will see below that (6.12) 



provides accurate results for incoming wave functions that lie outside of the scope of the 
theory in [T^, and where the transmitted wave function is clearly not Gaussian. 



7) A further advantage of (6.12) is that it only uses local information: the parameters Qc 



S and 7 are determined by the derivatives of the potential at the coupling point, and the 
incoming wave function is only needed at the crossing time. This immediately suggests 
an algorithm for computing non-adiabatic transitions, even for complicated, non-contant 
energy levels. As in surface hopping models, one evolves the wave packet on the initial 
adiabatic surface, until one detects a local extremum in the coupling functions 6'{qo{t)), 
where qo is the center of the wave packet at time t. One then com pute s the coupling 



functions k^^i from the local shape of the potential, and uses them in (|6.8l). The choice of 



the optimal n is of importance, and we suggest a way to find it in Definition [T] below. In 
(6.7), one can use the free propagator, since the transition region is small, and the energy 
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surfaces are approximately parallel to first order in q, as their distance is minimal at a 
crossing point. Then in a similar way as above, one gets an explicit transition formula. 



6.3. Numerics. We now show that (6.12) is in excellent agreement with highly precise 



numerical solutions of the full two-band Schrodinger equation. For solving the latter, we 
use standard methods, including a symmetric Strang splitting. We denote the projection 
of the numerical solution onto the lower eigenspace by (f)-{q, t). We compare this with the 
result given by (6.12), which we denote by V'- as before. We compare our results in the 
Fourier representation, calculating the inverse e-Fourier transform of (p- with a standard 
FFT. The final time t is chosen so that \\4>-{q, t)\\2 is constant under further time evolution 
in the exact calculation. 

For our potential we choose 

= f arctan(tanh(^f^)), p = S = 1/2 (6.14) 



This gives 9'{q) = c/2/ cosh(ag), with singularities closest to the real line at 
the residue at which gives 7 = — In particular, we take c = — vr/S, 
giving qc = 1 and 7 = 1/3. The choice of 6 over the case 9{q) = arctan(g/5) 



in ([L2 

ibig'c = ±i 
a = tt/2 

(which would make 9^ = in (5.2)) is made to increase the rate at which the potential 
becomes flat. This reduces the necessary computation time for the numerical solution. If 
anything, we would expect the asymptotic results for 9{q) = arctan(g/(5) to be better. 

Our first choice for the wavefunction in the intially occupied band is a Gaussian wave 
packet. At time t = 0, it is given by 

Vr;o'(p) = (2vr£)-i/4exp(-(P^). (6.15) 

As the crossing point is at x = by our choice of potential, ^^+,0 is sitting right at 
the middle of the crossing region. Since the eigenvalues are constant, it may be evolved 
backwards to to < exactly on the upper level: 

-1/4, 



^+ {p,to) = (27re)" 



'exp((^: 



exp{-Hoip'/2 + 5)). 



(6.16) 



We use (6.16) along with ■(/;_ (p, to) = as initial conditions for our numeric solution, and 



take to sufficiently negative in order for the wave packet to be well away from the crossing 
region. 



4 . 10-^ 

2 . 10"^ 

8 . 10--' 
4 . 10--' 
2 . 10--' 



10-' 
3 . 10-- 

10-' 
3 • 10"' 

10"' 
3 . 10-* 

10-' 
3 . 10-' 

10--' 



a) Po 



b) Po 



Figure 1. Relative errors between numerical results and (6.12), for a 



Gaussian wave packet, with different values of e and po', on a logarithmic 
scale. 



Our numerical studies show that (6.12) is in excellent agreement with the numerical 
solution for a wide range of e, ranging from as large as 1/10 to 1/50, at which point a 
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further reduction in e makes the numerically exact calculations very time-consuming. In 
Figu re [ij we show the relative error in the norms between the numerical calculation 



and (6.12 ), i.e. 



|2/||</>- 1|2- In each case, the step size in the numerical calculation 



was reduced until the difference between two subsequent numerical solutions was at least 
one order of magnitude smaller than the error to the solution obtained from (6.12). 

Figure [l] a) also shows that, as predicted in Section [5] ( |6.12 ) is not asymptotically 
correct for fixed p: After an initial increase in accuracy due to the decrease in e, the 
relative error becomes larger again as e decreases further. That this does not affect the 
practical usefulness of (6.12) becomes clear when we consider orders of magnitude: for 
Po = 2 and e = 1/50, we have HV'-lb ~ 6 * 10""'^'', which is much smaller than is useful in 
practice, while the relative error is still excellent at about 4 * 10~^, albeit deteriorating. 
On the other hand, for pQ = 2 and e = 1/10, we have HV'-lb ~ 0.014, which is certainly of 
a physically measurable size, and the relative error is still of the order 10~^. Finally, for 
e = 1/5 and po = 2, we have HV'-lb ~ 0.11, and the relative error is still below 0.03. We 
see that, as is often the case in asymptotic formulae, the actual error is much better than 
could be expected from the a priori error estimates: This is particularly obvious in Figure 
[l]b), where the relative error initially decreases like e~^^^ before saturating, while theory 
only predicts a ^/e decrease. Orders of magnitude in this case range from HV'-lb ~ 0.138 
for £ = 1/10, with relative error below 0.025, to HV'-lb ~ 6 * 10~^ for e = 1/50, with 
relative error below 2 * 10~^. 

Our second numerical example 
is an incoming wave packet that is 
strongly non-Gaussian. We choose 
^\p) = le-(P-Po)''/(^^) , where 

Z normalises the L norm of ■0+ • 
The potential V is precisely as in 
the other example, and we choose 
Po = 5 and e = 1 /50. Figure [2] 
shows the absolute value of the 
transmitted wave packet in the 
scattering region, in Fourier rep- 
resentation. The relative error in 
the L^-norm in this case is around 
5 * 10~^, similar to the one in 
the Gaussian case. In particular, 
the pointwise relative error is ex- 
tremely small in the region where 

- — - e 

tp-. is concentrated; thus, while 

Figure [2] shows the result of applying (6.12), the plot showing the numerical calcula- 
tion would be indistinguishable from the one given. Note also that the momentum shift 
relative to the energy conservation value \/pq~+~S^ ~ 5.02 is clearly visible. Obviously, 
the transmitted wave packet is strongly non-Gaussian, with a surplus of high momentum 
components having made the transition. 




Figure 2. Transmitted wave function 
for a non- Gaussian incoming 
packet. 



wave 



6.4. Optimal superadiabatic transition histories. From the time-adiabatic theory 
we know that in the optimal superadiabatic representation, the transmitted part of the 
wave function builds up monotonically, and has the shape of an error function as a function 
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of time. We are now going to show that a similar property holds for the Born-Oppenheimer 
transitions, after suitable modifications. 



Our starting point is Theorem (6.7), which gives a formula for transitions in the n-th 
superadiabatic representation. Our first task is to find a precise meaning to optimality of 
a superadiabatic representation. One natural idea would be to require that in the optimal 
superadiabatic representation the map i >— > / \ip-{x, t)pdx increases monotonously during 



the transition. But given the complexity of formulas (6.7) or (6.10), this condition is rather 



difficult to check. We will instead choose anoth er c ondition which, as we will argue, should 



be equivalent. The basic idea is that through (6.7), ip'L{k,t) is given by an integral where 
the integrand is expected to be both highly oscillatory and sharply peaked in s and k. 
Thus the main contribution occurs at points where the integrand has either a stationary 
phase, or a maximal absolute value. The locations of both of these depend on n, and if we 
want to have any chance of seeing a 'nice' transition history, we have to chose (if possible) 
n such that these locations coincide. Thus we define 



Definition 1. Let u = e*e^^ Jn+i^ ^/;^q, as given in (6.1). We write u{s,k,n) = 

exp{^{—X{s,k,n) + iY{s, k,n))) . Let (s*(n), A;*(n)) be the location of the minimum of X 
on the real line for given n. We say that n is an optimal superadiabatic representation if 
dsY{s^{n), k^{n)) 0. 

The ~ sign in the above definition means that in principle, n has to be integer, and thus 
equality might not hold anywhere. On the one hand, this problem should be less severe 
when e is extremely small, and on the other hand, usually Jn+i is given by an explicit 
formula, in which case we can interpolate and take n to be real. In that case, we will 
usually be able to fulfil dsY{s^:{n), k^,{n)) = exactly. 

That the above definition indeed gives 'optimal' transition histories can be seen by 
the following simple calculation. Let us write X^k = 5^X(A;*, s=k) etc., and k = k — k^, 



s^,. Expanding the exponent of (6.7) around k^ and s* then gives 








Thus for k = 0, i.e. at the maximum of the transmitted wave function, we have 

i „ _-^^n li-X(k,,s,)+iY{k„s,)) ( ( /X _IY 

ei*^- AK^) ^ 1 ,^ 1 + erf (t - s.). ^^'^ ' 



2^Xss - \Yss \ \ V 2e J J' 

We see that the transmitted wave function, when adjusted for the propagation in the lower 
band, has the shape of an error function at its maximum /c*. The only unusual feature 
is that this error function is actually evaluated along a 'diagonal' in the complex plane, 
rather than on the real line. Nevertheless, the resulting shape will be close to monotone 
unless Yss is much larger than Xgs, in which case there are some oscillations around t = 0. 
We do not know whether this case is likely to happen in practice. For k ^ 0, the error 
function behaviour deteriorates, but since V'-,n is a semiclassical wave function, we are 
only interested in ~ ^/e, in which case the behaviour is still similar to the one at the 
maximum. 

We now want to show that optimal superadiabatic representations exist in particular 
cases. As an example, we pick again the situation of constant eigenvalues. We rewrite the 
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integral in (6.10) as 



ds / dr/exp 



-M{k, 7]) + i(j){r]) + -s{k^ -rf - 4(5)) 



(6.17) 



Here, e~^^^^ is the combined modulus of the transition kernel originating from J^j^i and 
the wave function, and (f) is the phase of the wave function, hence depending only on r/. 
We can treat more general forms of the coupling function k~ (p, q) than the one leading 
to (6.10), including the full coupling function with all powers of p included. The only 



requirement is that, to leading order, k should be symmetric, so that its Fourier transform 
is real. This is generically true: high derivatives of the pair of first order poles in the complex 
plane determines the shape of k, and these are either symmetric or antisymmetric, giving 
either purely real or purely imaginary Fourier transforms. 

Given (6.17), we now let r/* = r/*(n) and A;^, = k^:{n) be the place where M{k,r]) is 
minimal. We expand M to second order around (/c*, ry*). In order for the phase to be also 
stationary, we need dr^(l){rj^:) = srj^, which is the equation determining the transition time 
We should keep in mind that any further explicit calculations will only make sense for 
s close to Now, we expand the phase to second order around r/* as well, and compute 
the resulting Gaussian integral in rj. The result is 



;g-^(M(r,*A)-i</'('?.) + |A4fcP) ^ 



27re 



2^ , is{fc2+2fctfc+fc2_^2_45)+. 



-{l»{s-s*)- 



(6.18) 



ds. 



Above, we use again the notation k = k — k^,, Mj^f^ = d^M{k^,,r]^) etc. As before, we 
concentrate on the case A; = 0. Let us assume that we can find n such that kl — rjl — Ad = 0. 
It is then easy to check that for the remaining integrand, both real and imaginary part of 
the exponent are stationary at s = s*, and that s* is a maximum of the real part. Thus 
for constant eigenvalues, we can find an optimal superadiabatic representation if we can 
solve the equations 



dkM{k,r]) = ^, d^M{k,ri) = 0, k^ - rf = A5 



(6.19) 



simultaneously, and if the resulting pair {k,rj) is a minimum of M. Note that as the above 
equations also depend on n, we are solving a system of three equations with three free 
variables, which means that we can hope for a solution. Note also that the third equation 
is connected to energy conservation: r/ is the incoming, k the outgoing momentum and 5 
is the gap between the energy levels. 



We specialise further in order to solve (6.19). We choose the upper band wave function 



to be a Gaussian wave packet with momentum p^, so that 



1 



(cj2^e)i/4 



exp 



1 



iv-Po)' 



(6.20) 



Choosing ip+fi {rj) real-valued amounts to choosing the packet to be at x = at time 
t = 0; putting the avoided crossing at x = in addition ensures that the transition occurs 
at time s* 



0. Now, using the integrand in (6.10), we get, to leading order. 



M{k, T]) = -ne(ln(r - r]^) - ln(4(5)) + qc\k - ??| + 



2cj2 
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By the third equation in (6.19), we know k > r], which removes the absolute value above. 
Taking derivatives, 

k 



= 2ne- 



rj 



1^2 _ ^2 ' 

which together with the equation k'^ 
k = \Jrf- +4(^, 



+ 







-2nE 



k'^ 



+ 



k 1 



4(5 lead to 
V-Po 



n 



26qc 

ek 



(6.21) 



The first two equations are independent of n and e, and determine rj^,k.t. This is a spe- 
cial feature of the constant eigenvalue situation, and should not be expected in general. 
The third equation determines the optimal superadiabatic representation. Note that n is 
connected to the optimal superadiabatic n for the time-adiabatic situation: there, we have 
nta = 25qc/e, with the 'momentum' (i.e. speed on the time axis) normalized to one. But 
the tricky point that can't be easily guessed is which value of k to pick in the formula for 
n^,] the naive guess of using the incoming momentum po would be totally wrong. The more 
sophisticated guess of using the mean momentum at the crossing point, (?]* + A;=k)/2 would 
be closer for small (5, since the true value fulfils k^, = ^{k^ + ^/rf + 4(5); but it would still 
be far off for finite (5, which are of main interest here. 





ll^-.0(-,t)l|2 



IIV'-,l(.,t)l|2 



IIV'-.2(.,t)||2 





ll^-,4(.,t)l|2 



liV'-.5(.,i)|l2 



Figure 3 . Time development of the norm of the transmitted wave func- 
tion, from the adiabatic up to the fifth superadiabatic representation. The 
dashed line is the theoretical prediction of the optimal superadiabatic tran- 
sition history. 



We close by comparing the effective formula (6.18) to numerical results for the transition 



histories in various superadiabtic representations. In order to obtain the optimal supera- 
diabatic representation at both low and close to integer values of n, we have to choose 
the parameters somewhat carefully. Our choices are the potential (6.14) with c = — 7r/3, 
27r/5 and 5 = 3/32. This gives qc = 5/4 and 7 = 5/12. We take e = 0.02923, and 



a 
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in the incoming wave function (6.20), we take o"^ = 2 and po = 2.5. Solving (6.21) then 
yields t/* 2.57, k^: « 2.64, and « 3.04. Thus the optimal superadiabatic representa- 
tion should be the third, and this is clearly confirmed by numerical simulations. Figure |3] 
shows the norm of the transmitted wave function (calculated by a numerical solution 
of the Schrodinger equation), in the adiabatic and all superadiabatic representations up to 



n = 5, as a function of time t. The dashed line is the prediction of formula (6.18) with ry* 
and fc* from above inserted. We see that indeed the optimal superadiabatic representation 
is at n = 3, while below n = 1 and above n = 4 oscillations grow. The reader should 
note the similarity with the plots shown in |22] . This is another confirmation that the 
time-adiabatic approximation is appropriate to qualitatively understand the mechanism 
of non-adiabatic transitions, as far as the population of the lower level is concerned. On 
the other hand, in order to obtain quantitatively correct results, or more detailed informa- 
tion about the transmitted wave packet (e.g. momentum spread of phase shift), the full 
quantum mechanical treatment, as given in the present work, is indispensable. 
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